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CN . Abstract 

The derivation of the Navier-Stokes equation starting from the Liouville equation using projector 



o 



techniques yields a friction term which is nonlinear in the velocity. As has been explained in the 1. 
version of this paper, when the second-order part of the term is non-zero, this leads to an incorrect 
formula for the equation. 

In this 2. version, it is shown that the problem is due to an inadequate treatment of the correlation 
function K%. Repeating the calculation leads to zero second-order part. The Navier-Stokes equation 
i-C ■ is correctly derived by projection operator technique. 
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<D 1 Introduction 

bp: 

The derivation of hydrodynamic equations by Zwanzig-Mori technique is a method well-established in 
the literature. In version 1 of this paper, the author reported a problem which occurred when the Navier- 
r/3 ■ Stokes order of the momentum equation is considered: As is well known, the non-dissipative part shows 

to be of second order in the fluid velocity, u. The dissipative part is also non-linear in u. In order to 
obtain the correct form of the Navier-Stokes equation, it is necessary for the second-order part of the 
dissipation term to vanish. As far as the author knows, this aspect has not been investigated earlier. The 
calculation in version 1 yielded a second-order part different from zero. 

In the present version it is shown that the problem resulted from an incomplete formula for the 
■ correlation function K2- Actually, the second-order term is zero. Thus, the Navier-Stokes equation is 

\ correctly derived by projection operator technique. 
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O ' 2 Hydrodynamic equations 
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For the basic definitions of microscopic variables, readers are referred to [T]. In order to keep the paper 
reasonably self-contained, the results of the Zwanzig-Mori analysis are copied from there. The analysis is 
for incompressible constant density/temperature processes. The hydrodynamic momentum equation is a 
■ specialization of the general mean value equation (3.1): 

p{-£ + u c V c u a ) = -V a P + D a (2.1) 
at 

In (3.15) this equation is given for stationary processes. Latin indices run from 1 to 3. p is the fluid 
density; u as well as the pressure P and the dissipative force D depend on space and time. For the latter 
the formula is obtained: 

D a (x, t) = /3V C J dt 1 J dx'i? abc(i (x, t, x', t')V> 6 (x', t') (2.2) 
with the kernel function: 

i? abcd (x, t, x', t') = ([G{t', i)s ac (x, i)]s M (x', t')) L ,t> (2.3) 
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(3 is the inverse kinetic temperature. () L,t denotes the expectation with respect to local equilibrium, 
cf. pQ (3.4) to (3.6). s ac is the projected momentum flux density: 

Bacfct) = (1 - P(*))s ac (x) (2.4) 

s ac is the momentum flux density. V(t) is the Zwanzig-Mori projection operator; for any phase space 
function g, it is defined: 

Vg = {g) L + (gSa) L * (Sa6a)^ 1 *8a (2.5) 

For shortness, the time parameters have been omitted here. * denotes an operation which consists 
of a product, a summation over 5 index values and a space integration, a are the microscopic densities 
of the conserved quantities. Sa = a — (ci)l. denotes the inverse matrix. Q(t',t) is a time-ordered 

exponential operator: 

G(t',t) = exp_{J dt"£(l - V{t"))} (2.6) 
C is the Liouville operator, cf. pQ (2.4). 



3 2 nd order term of the friction force 

In this paper, we consider the second-order verlocity approximation of the friction Term D. This should 
not be confused with the second-order wave number approximation of the linear part of D which leads 
to the well-known Stokes form of that term. This latter approximation is not considered here. 

As is explained in pQ, the 2 nd order part D^ 2 ) of the dissipative force is obtained from (|2.2j) by inserting 
the linear part i?*- 1 ** of the kernel function given by p] (4.1), (4.2). These formulas have to be generalized 
for time-dependent processes: 

Kbcd( x > x > t , t ) = -P J o dt J dx t „) |u=oMx ,t ) (3.1) 

b e = ~(3u e is the momentum part of the conjugated parameters in the formula for the local equilibrium 
probability density, see [1] (3.6). The functional derivative is calculated in the appendix. 
The second-order term of the friction force consists of a parametric and a functional part: 

D 2 = D 2p + D 2f (3.2) 

When the parametric part of the derivative (|A.7[) is inserted into (|3.ip . and the result ist introduced 
into (|2.2[) . one obtains the parametric part of the second-order friction force. We find: 



(D 2p ) a (x,t)=[3 2 V c J dx' J dx"fdt' ([e( 1 -^) £ (*- t ')(s )ac(x)](so) M (x')p e (x"))ou e (x",t')V> b (x',i') 

(3-3) 

()o denotes total equilibrium expectation; So is the flux variable projected with (1 — Vq); Vq is the 
total equilibrium counterpart of V (|2.5p . - If one inserts the A th part of the derivative (|A.16[) into (|3.1|) , 
one obtains the functional part of the friction force: 



(D 2f ) a {x,t) = (3 2 V C J dx! J dx" dt' £ dt"x 



x ({^e(i-^Wt''-t') PoPe{x »)e^oW~^ (3.4) 

It is necessary to derive formulas for the appearing correlation functions in order to get a definite 
result. - We begin with transferring the formulas into Fourier space. In addition, within the correlation 
formula, we switch to orthonormal variables h e , r ac , r ac (see [T] (4.7)) . 

The hydrodynamic velocity varies over much larger time intervals than the correlation functions. 
Therefore, a Markovian approximation of the process is applied which results in the approximation for 
the kernel function f(t) in the integrals: f(t) —> 5(t) J Q dt' f(t') . For the parametric part, we obtain: 

(D 2p ) 1 (t)=^p i j^Te J d( lJ dq'(N p ) 123 u 2 (t)u 3 (t) (3.5) 
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(N P ) X23 = dt'de^-^'i^f^qfr^o (3.6) 
Jo 

Number indices are introduced which have been used by many authors. The number is a combination 
of an index and a wave number. Certain numbers are reserved for a specific wave number, as is shown in 
the following table: 

k: 1, 4, 7,.-. 
q: 2,5,8, ••■ 
q' : 3,6,9,. 

Number indices appearing pairwisc include a summation over the original indices; on the other hand, 
wave number integration will always be shown explicitly. In (|3.6|) . (ikf)i is written for ik^fid- - N p is 
identical with the quantity N in [T] (4.10), which is calculated there. The result (4.27) is the "reduced" 
form of the kernel function which means that the factor (2-7r) 3 (5(k — q — q') is excluded. The full formula 
reads 

(AW 123 - (2vr) 3 (5(k q q)~ik d S 123d (3.7) 
Details of the matrix S are not needed and are therefore not repeated here. - (|3.4p changes to: 



(D2f)l(t) = 




The projection Vq in the correlation is performed: 



(U 2 /)l(t) = ~P h P i J^ J d <lJ d <l' j dt> j t , dt " 



x([e( 1 -^) £ ( t - t ''H*kr) 1 ]^^) ([| 7 e( 1 -^) c ( t ''- t ')/ l2 ](jqf)*) oU5 (^)«3(t'') (3.9) 

() denotes the 3-point-correlation without the factor (27r) 3 <5(k q q'), which we call the reduced 
correlation function. - Finally, symbols are introduced for the correlations appearing in 



(D 2f ) 1 (t) = -(3ipi-^ /dq/dq' J'dt" J* dt\K 3 j 123 (t-t") diK2)25 Jt" f) ' u 3 (t")u 5 (t') (3.10) 

(K 3 ) 123 (t) = ({e^o)ct {{kf)i]h * h * )o (3 n) 

(K 2 ) 25 (t) = <[ed-W h 2 ](iqr)S)o (3.12) 

4 Correlation function K<i 

We start the inverstigation of K 2 by calculating its initial value: 

(#2)25(0) = (h2(iqf)* 5 )° = ([(1 - V)h 2 ]{i<tf)t) Q = (4.1) 

Here we have used the general properties (/(l — P)g)o — ([(1 — T'O/K 1 — ^ > )5)o and (1 — V)h = . 
Next we derive a relation between K 2 and correlation functions which are defined with the non-projected 
exponential operator e ct . For K 2 , we use the relation: 

(ikf)i = —hi —iuJuhti (4-2) 

We obtain the decomposition: 
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(K 2 ) 25 (t) = -([e( 1 -^) £t / l2 ]^)o+*c.58([e (1 - T,o)£t / l2 ]^}o 

-(#21)25 (t) + iu> 5S (K 22 ) 2S {t) (4.3) 

w 25 = ((cf)2^5)o: and the symbols in the second row are denotations for the quantities in the first. 
We use the operator identity |T] (4.11): 

c (l-V a )Ct = e Ct _ f* dt , e Ct' Vq£ e (l-p )£(t-i') (4 4) 

JO 

Application to K 2 i yields: 



(*ai)«(t) = -/ dt'{\e c >"v„£< ! ' 1 - TM(, - , ' ) h 2 ]hl)„ + {\>!"h 2 ]ht)„ 

= -Ki^-tm^-^^. (4.5) 



dt' dt 

This is an integral equation for K 2 \ in terms of the non-projected duple correlation ^2)25 (£) = 
([e Ct h 2 ]hl)o- (C 2 ) means (C 2 ) without the factor (27r) 3 <5(k — q), which we call the reduced form of the 
two-point correlation. C 2 obeys the equation [T] (4.24): 

d{C 2 ) 25 (t) . . 

J t = -K28(C<2)B5W > t>0 ( 4 - 6 ) 

K28 = «W25 +7 25 (4.7) 

7 is the dissipation matrix, the space-time integral over the memory function of the process which 
is defined with the linear projection operator employed here. In .2] it is emphasized in connection 
with formula (3.18) there, that in general one has to distinguish this from the corresponding quantities 
defined with the multilinear projection operator defined there. For the limited purpose of the present 
paper, this difference can be ignored. The condition t > in (|4.6[) is essential, since at t = the time 
derivative of C 2 is discontinuous: While from the definition of C 2 and the conservation relations we find 

■^-gt \t=o = —ioJ 2 5, from (14. 6p one concludes d ^ 2 ^ 25 |t=o + = — ^25- - (|4.6| is introduced into (|4.5|) . and 
the integration parameter t' is transformed: 

(K 21 ) 25 (t) =Ks,u [ d*'(^2l)28(0(C2) 11 5(i-0+«2,ll(C2)ll, 5 (*) , *>0 (4.8) 



Differentiation with respect to time and insertion of (I4.8l) . (|4.6p yields: 

d( ^ 25(t) = -Ku,5 {[(K 2 i) 2M (t) ~ «a,ii(Ca)n,u(*)] + ^11(^2)11,14 (*)} + Kss^aiM*) 

= , t > (4.9) 



Thus, for t > 0, i^2i is a constant. The value ist found from ([4.8)1 by taking the time limit t = 0+: 

(#21)25 = K25 , <>0 (4.10) 

For calculating the properties of i^22, we need an operator identity similar to ()4.4[) : 



(i-v)ct = e ct_ dt , e {i~v)ct' VCe c{t-t') (4n) 



e v 

Both formulas can be found from more general identities in [3]. The calculation ist quite similar to 
that of K 2 i, and the result is: 

<^M) =0; t>0 (4.12) 
(#22)25 = S 25 , t>0 (4.13) 
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Introducing these results into (|4.3p , we obtain: 



^t 5(f) =0. ^>0 (4.14) 
at 

(K 2 ) 25 = - 725 , t > (4.15) 



Combining (14. 15)) and (|4.1j) yields: 



(^2)25(t) = -725©W . *>0 (4.16) 
O(i) is the step function. By Differentiation: 

^ = -7 2 ,^W, t>0 (4.17) 

It is understood that the dimension of the delta function corresponds to that of its argument; always 
the same symbol S is used. The factor 2 is necessary since the delta function has its peak at the boundary 
of the definition range of K 2 . - Formula (|4.17[) constitutes the essential difference to the first version 
of this paper, where the result (4.9) stated d(if2 j t 25(t) = , t > 0. By inserting (|4~T7|) into l[3~TUj) and 
evaluating the delta function: 

(D 3f )i(t) = y rfq y dq' J dt'(K 3 ) 153 (t-t')u 2 (t')u 3 (t') (4.18) 

Since the hydrodynamic velocity varies on a much longer time scale than the equilibrium correlation 
function, it is again possible to apply a Markovian approximation: 

(£>a/)i(t) = /3 V j d *j dq'(N f ) 123 u 2 (t)u 3 (t) (4.19) 

/•OG 

(JV»i23 = -7 25 / dt'(K 3 ) 153 (f) (4.20) 
Jo 

5 Calculation of 



In contrast to the first version of this paper, we now need to calculate the correlation function K 3 . We 
again use the identity (|4.4p ; with a calculation which parallels pQ (4.10) to (4.12) we obtain a formula for 
K 3 : 

(K 3 ) 123 (t) = -7i4<[c £t h A ]h* 2 h* 3 ) + ([e £ *(zkf ) 1 }h* 2 h* 3 ) Q (5.1) 
The first term contains the non-projected single-time triple correlation function: 

(C 3 ) 123 (i) = ([e ct hx]hlh* 3 ) (5.2) 
The second term can be transformed by the relation (|4.2p . We obtain from (|5.ip : 

(K 3 ) 123 (t) = -«u(C 3 ) 423 W " d{C ^ [t) (5.3) 

We switched to reduced correlations since this is more suitable for the following calculation. - As has 
been stated in pQ, for C 3 expression (4.20) has been derived by multilinear mode coupling theory: 

(C 3 ) 123 (t) = (C 2 ) 14 (t) J 423 - / dt'(C 2 ) 14 (t - tyk d S 45 UC2) 52 (t')(C 2 ) 63 (t') (5.4) 



JO 

J is given in pQ (4.23). For K 3 , we then find: 

(K 3 ) 123 (t) = ife d Si 5 6 d (C 2 ) 52 (i)(C 2 ) 63 (t) (5.5) 
The solution of (|4.6p to the initial condition (C 2 ) 14 (0) = S14 is: 
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(C a ) u (t) -e"^ 4 (5.6) 
This is inserted into (|5.5p ; by integration we obtain for Nf (|4.20|) : 

(-^/)l23 = -725(^58^36 + S 58 K36)~ 1 ikdSl86d (5.7) 

For insertion into (|4. 19[) . we may replace (Nf) 123 by the symmetric form ^((Nf) 123 + (Nf) 132) which 
we denote by the same symbol: 

( N /)l23 = ~ 2(l25^G + ^25736) ( K 58<^69 + ^58'«69)~ 1 «fcrfS'l89rf 

= --jikdS\23d + ~(iw 2 5#36 + <525«W 36 )(/C 58 (5 6 9 + ^58K69)^ 1 «fc ci S'i 8 9 ci (5.8) 

When this is introduced into (|4.19p . it is seen that the second term vanishes because W25 = ^5 92 (&5 
certain constants), and the incompressibility condition reads 92^2 = 0. Thus, for the relevant part of Nf, 
we have: 

(JV»i23 = -(2^) 3 <5(k q q')^ik d S 123d (5.9) 



When we compare with (|3.7|) . we obtain from 

D 2 =0 (5.10) 

Thus, with the recalculation of K2 (|4. 16[) . we now obtain the result that the second-order part of the 
friction force vanishes. 

6 Stationary processes 

We shortly turn to the earlier paper pQ, where we considered stationary processes. When the calculation 
is performed in Fourier space and orthonormal variables are introduced, inserting (A9) into (4.1) and 
importing this into (3.2) of that paper yields the consecutive formula fur -D2/: 

(D 2f ) 1= pip§ JdqJctfk^J & J dt''([—e( 1 -^ a 'vh^ 1 -^ a ''hc]r* 2d ) u 2 u 3 (6.1) 
As in 1], the projection operation is performed: 

(S 2 /)i=/3M^ /dq/ ^ ^ rft "<[ eC1_7 ' >jCt,/ ^id^S^Do J7<[ e<1 ~ 7 ' )JCt, ^6]^d>o« a «3 (6.2) 

It has been taken into account that we have: 

{[c (i-v)ct" hc]h * )o = {[e c(i-v)t" hc]{1 _ v)hl)Q = (6 3) 

Definitions (|3~TTj) . (|3~T2|) of the present paper are introduced: 



1 r r r°° r c 

(D 2f )i = -/3M 7^)6 J dq J dq'u2U 3 J dt" (K 3 ) 153 (t") J 



dt 



,d(K 2 ) 52 (t>) 



52 V 
dt" 



= ~^ 2pi J2^ J dq J dq ' U2U3 L dt '^ K ^ 3 ^[}^J K ^(t')- (^2)52(0)] (6.4) 

It has been assumed in pQ that for large t the factor variables of if 2 become statistically independent; 
this led to lim^^oo (^2)52(0 = 0. But as we see from (|4.15|) . for K 2 this assumption is not true. When 
the correct formula is incorporated, the calculation leads to the final result (|5.10D . 
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7 Summary 



As is stated in the introduction, D2 = ensures that the Navier-Stokes equation is obtained correctly 
as second-order velocity approximation of the mean momentum equation derived by projection operator 
technique. 



A Appendix: Calculation of the functional derivative 



The kernel function (12. 3|) depends on b e 4-fold, namely, in the formula for the local equilibrium density 
p L , inV contained in s ac and §bd, and in the operator Q. For abbreviation, the formula for the derivative 
is written: 



SR abcd (x,x',t,t') 
Sb e (x",t") 



E 



SR 
5b 



i(0 



(A.l) 



In the first three of these, the time dependence of R is parametric via the time dependence of p L , 
at the time instant t or t' . In these cases, the time-dependent derivative is equal to the corresponding 
stationary derivative times S(t" — t) or 8{t" — t'). The stationary derivatives are given in [1] (A3) to 
(A5). We still have to take these formulas at u = 0; then, p L switches to p , the probability density 
of total equilibrium (for a given density and temperature), and therefore Qlj to ()o , the equilibrium 
expectation; V reduces to Vq , the corresponding total equilibrium projection operator; s ac to (so)ac, the 
flux density projected by Vq\ Sp e to p e , and G(t', t) to e^-f^cit-t )_ -y^ e obtain: 



SRabcd 


(1) 


Sb e 


u=0 


SRabcd 


(2) 


Sb e \ 


u=0 


SRabcd 


(3) 


5b e J 


u=0 



We have: 



= -S(t" -t')/3([c £(1 - T ' o)( *- t ' ) (so) QC (x)K(x")(so) M (x'))o 
= 6(t" - t')(3 ([e £ ( 1 ^)(*- t ')(s ) QC (x)]P p e (x")(so) M (x'))o 
5(t" - t)P ([e^-'PoKt-t') ^ oPe ( x ")(j Q ) ac ( x )](s ) bd (x'))o 



{[e C(l-Vo)t VoA]{1 _ Vo)B)o = ( ^ 1 _ Vo)e CH-Vo)t VoA ^ 1 _ Vo)B)o 

= 



(A.2) 
(A.3) 
(A.4) 



From this we find: 



SRa 



bed 



(3) 



u=0 



(A.5) 



(A.6) 



The remaining two parts yield the parametric part of the functional derivative: 



SRabcd 




SRabcd 


(i) 


SRabcd 


(2) 


[ 6b e 


V 


[ 5b e 


+ 

u=0 


[ 5b e 


u=0 



-*(*" - 0/3 ([e £(1 ~ 7 ' o)(t - t ' ) (so)ac(x)](l - 7> )p e (x")(soMx'))o 
-6(t" - 1')(3 ([e( 1 - T '«) £ ( t -*')(s ) ac (x)](s ) M (x')p e (x"))o 



(A.7) 



Note that in the 2 nd step, in the exponent the sequence of operators is reversed; compare (|A.5[) for a 
similar operation. 

The 4 th part of the derivative is defined: 



SRa 



bed 



6b e (x.",t") 



(4) 



5Q{t',t) 



6b e (x", 



frAi - P(i)Kc(x)](i - V{t'))s bd {x!))L,« 



(A., 



By temporarily approximating the Integral in (|2.6p by a sum, one finds for the derivative of G(t',t): 
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sg(t',t) 



8b e (x.",t") J t , " v ' ; 56 e (x",i") 
Moreover, for the middle factor of the integrand which depends parametrically on r, we have: 

S(C(l-V(r))) 



b e (K",t") 



8{t" - T)£V(t")5 Pe (x",t")(l - V(t")) 



This is introduced into (|A.9[) . The following structure arises: 



dTS(t" ~ T)F(t",r) = 0(i" - t')0(t - t")F(t", t") 



with <d(t) being the step function: 



We obtain: 



5b e (x",t") 
This is introduced into (|A.8|) 



(A.9) 



(A.10) 



(A.ll) 



(A.12) 



e(t" - 0®(* - t n )g{t' ,t n )cv{t n )5p e {^' - v{t"))Q{t" ,t) (A.13) 



60 <7 



(4) 



For u = 0, this reads 



= [3Q(t" -t')Q{t-t")x 

x ([£?(*', t")CV{t")8p e {x" , t")(l - V(t"))g(t", t)&ac(x, t)]hd(x', t')) L , t , (A.14) 



SRabcd 



5b e (x",t") 



(4) 



u=0 



x ([e^-W-O £7> oPe (x")(l - 7>„) e £ ( 1 - T ' ^ t - t "))(so) ac (x)](so) &(i (x'))o (A.15) 
After some manipulations, the result is: 



5Ra 



bed 



5b e (x",t>' 



(4) 



u=0 



= p@(t" -t')Q(t-t")x 

x ([ e (l-^) £ ( t "-*')(i _ T'o)^oP e (x")e( 1 - T ' )^ t -*"))(so) ac (x)](^) M (x'))o 

= -/3e(t"-t')e(t-t")x 

x([J 7 e( 1 -^)^''-Op Pe(x'Oe( 1 -^^( t -*''))(so)ac(x)](so) 6d (x'))o (A.16) 
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